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We introduce a probability current analysis of excitation energy transfer between states of an open 
quantum system. Expressing the energy transfer through currents of excitation probability between 
the states in a site representation enables us to gain key insights into the energy transfer dynamics. 

It allows to, i) identify the pathways of energy transport in large networks of sites and to quantify 
their relative weights, ii) quantify the respective contributions of unitary dynamics, dephasing, and 
relaxation/dissipation processes to the energy transfer, and iii) quantify the contribution of coherence 
to the energy transfer. Our analysis is general and can be applied to a broad range of open quantum 
system descriptions (with coupling to non-Markovian environments) in a straightforward manner. 


I. INTRODUCTION 

In the description of electronic excitation energy transport in systems of coupled sites - such as molecular aggre¬ 
gates, biological and artificial light-harvesting systems, coupled quantum dots - knowledge about the contributions 
of the different processes involved is important both for understanding and modeling of the energy transport mech¬ 
anisms in nature, as well as for design of artificial and bio-mimetic systems that enable efficient transport. The key 
processes involved in the energy transfer that are essential for this analysis are unitary dynamics, dephasing, and 
relaxation/dissipation. It has been found that the right balance between these processes can enable highly efficient, 
directed energy transport m- Another key aspect in the investigation of design-function relationships in natural and 
artificial light-harvesting systems is the identification of the specific pathways of energy transport in large networks of 
sites and the quantification of their relative weights 0H2- In such systems and for open quantum systems in general, 
it has long been appreciated that coherent dynamics play a role in the overall energy transfer, and the coherent 
features of transport have been intensively studied in a broad range of open quantum systems (10HT8j . In recent 
years, experimental evidence for long-lasting coherence in biological photosynthetic systems obtained from ultra-fast 
spectroscopies has also sparked questions and discussions about the contribution of coherence to the energy transport 
that drives photosynthesis [ 16] 'T8l - i28| . 

To address these aspects of energy transfer in natural and artificial systems, we introduce here an analysis in which 
the energy transfer between the states of an open quantum system is expressed through excitation probability currents. 
In this paper we undertake our analysis in a site basis that consists of states in which the excitation is localized on 
single sites, since we are interested in applying the probability current analysis to spatial energy transport. However, 
the excitation probability current analysis is general and is applicable to description of the energy transport in any 
basis. 

Describing excitation or charge transport by means of probability currents is commonly encountered in a number 
of situations in condensed matter [29l - iJT] and in electron transfer reactions [32]. However, such a description is 
usually not applied in the study of light-harvesting systems, where energy transfer has instead usually been evaluated 
in terms of time-dependent populations of excited states (see e.g. Refs. 0 Hi m E3H3H])- In our analysis, we 
obtain the probability currents from a continuity equation for the excitation probability of the sites. Since the total 
excitation probability inside the (open) system is usually conserved, such a treatment based on continuity equations 
and probability currents is widely applicable. In this paper we show that this analysis enables us to address the 
following tasks in a straightforward manner: 

1. Calculate the excitation probability currents between the individual sites, i.e., their direction and magnitude, 
thus revealing the pathways of transport. This can be performed for a large, complicated network of sites, where 
the excitation of the sites couples to either a Markovian or a non-Markovian environment. 

2. Quantify the respective contributions of unitary dynamics, dephasing, and relaxation/dissipation to the currents. 

3. Quantify the contribution of coherence (versus the contribution of populations) to the currents, i.e., to the 
energy transfer. 
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Here we shall understand coherence to be quantified by the off-diagonal elements of the system’s density matrix in 
the basis under consideration. As noted above, in this work we employ the site basis in order to study the spatial 
transport of energy. 

Our treatment is applicable to the tight-binding and Holstein-like models - whether described as open or closed 
quantum systems - that are typically used for modeling excitation transfer systems of the form of biological or artificial 
light-harvesting systems [39ll44j . In the probability current framework, we obtain formulas for the contribution of 
the different processes - coherence, dephasing, relaxation - to the currents (i.e., to the excitation energy transport) 
from the evolution equation for the system’s density matrix. Inserting the elements of the time-dependent system 
density matrix into these formulas, where the matrix elements are obtained either from a numerical simulation or a 
tomographic reconstruction from experimental measurements, one can then numerically evaluate the contributions to 
the currents for specific model situations. The analysis developed in this paper is applied in a companion paper [4] to 
numerical simulations investigating excitation energy transfer and design-function relationships in the Photosystem II 
protein complex that drives photosynthesis in higher plants. 

In the following, we first briefly review the description of energy transport in an open quantum system by means 
of a non-Markovian quantum master equation (Section [TT|) to provide the underlying basis for our probability current 
analysis. Then, in the main part of the paper in Section m we introduce our probability current analysis, where we 
express the energy transport through probability currents and derive the different contributions of unitary evolution, 
dephasing, and relaxation to the currents based on the quantum master equation. Using this framework, we then 
explicitly analyze the contribution of coherence to the energy transport. To illustrate the interplay between the 
probability current and the coherence, we analyze the dynamical equations for the case of a Markovian environment, 
for which the equations are particularly simple and yield useful general insights. A quantitative analysis for a (33 
site) system with a more general non-Markovian environment is given in Ref. j4j. Concluding remarks are given in 
Section [TV] 


II. QUANTUM MASTER EQUATION DESCRIPTION OF ENERGY TRANSFER 

Typically, in light-harvesting systems and other systems of coupled sites featuring electronic excitation energy 
transfer, the electronic excitation also couples to other degrees of freedom. For example, in molecular systems such as 
pigment-protein complexes in biological light-harvesting apparatuses, electronic excitation of the molecules (pigments) 
couples to internal vibrational modes of the molecules as well as to vibrational modes of the protein scaffold in which the 
pigments are embedded [431 i46j . This coupling to vibrations is often crucial in the design-function relationship, since it 
induces dephasing and vibrational relaxation which help to efficiently direct the energy transport to a target location 
(in biological systems, this might be a reaction center) where the energy is trapped and transformed pi l4l l37l 147] . It 
is therefore important to include such vibrational degrees of freedom in the modeling and simulation of the excitation 
transport, in order to take these important dephasing and relaxation effects into account. However, this can easily 
result in a large number of interacting degrees of freedom, making it challenging if not impossible to numerically 
treat all of these degrees of freedom explicitly. It is thus very common to use an open quantum system description 
in which the problem is divided into three components. These are, i) the “system” that contains the degrees of 
freedom of interest, i.e., for electronic excitation energy transport usually the excited electronic states of the sites, 
ii) the “environment” - usually all the vibrational degrees of freedom - and iii) the interaction between system and 
environment, i.e. usually the electron-vibration coupling [4211481 . In this open quantum system approach one then 
solves an effective evolution equation in the small space of the system degrees of freedom, where the system degrees 
of freedom are treated explicitly and the environment degrees of freedom are taken into account only implicitly. This 
approach makes even large networks of sites numerically manageable [38l [49]. One widely used class of effective 
evolution equations derived from an open quantum system approach is represented by quantum master equations. 
These describe the time evolution of a reduced density matrix of the system, which formally corresponds to tracing out 
the environment degrees of freedom in the density matrix of the total problem EMEUS!]. Thus, for the simulation of 
electronic energy transfer, the reduced density matrix usually contains the electronic degrees of freedom, i.e., where 
formally the partial trace over the vibrational degrees of freedom was carried out ®mm, or alternatively, the 
electronic and selected vibrational degrees of freedom [52] . 

Since such quantum master equation approaches are widely used for the description of energy transfer and are 
particularly popular for the modeling of light-harvesting systems [3 0 HE EU Eg, in this paper we shall use a 
quantum master equation description as the basis for our probability current analysis of the energy transfer. In 
this section, we therefore briefly outline the quantum master equation description, before introducing the probability 
current analysis in Section [HI| 

We consider electronic excitation energy transport in a system of coupled sites, and as a dynamical model consider 
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a quantum master equation of the form 

d t p(t) = V(p(t)) = —i[H, p{t)} + C non—unitary (p(^) ); (1) 

that describes the time evolution of the reduced density matrix p{t) of the system, which contains the electronic 
degrees of freedom of the sites, with initial condition poj with H the system Hamiltonian that contains the energies 
of the sites and the couplings between these. (Here and throughout the paper we set h = 1). The first term in 
Equation ([l]) describes unitary dynamics of the system and the second term describes the non-unitary contribution 
to the system dynamics due to interaction of the system with the environment. The probability to find electronic 
excitation localized on site n at time t is then given by 

Pnn(t) = (n\p(t)\n), (2) 

where \n) is the state in which only site n is excited and all other sites are in the ground state. Thus, the states | n) 
span the one-excitation subspace. The transport of electronic excitation energy between the sites will manifest itself 
in the time-dependent behavior of these excitation probabilities p n n{t) of the sites. 

We assume that the non-unitary term of Equation ([T]) is a sum of two contributions, one of dephasing and one of 
dissipation/relaxation between electronic states due to the coupling to the environment: 

£non-unitary(p(t)) = £° ephaS (p{t) ) + C^ipit)) . (3) 

We next consider the form of these dephasing and electronic relaxation contributions. For our dynamical model, we 
shall assume that Equation ([TJ is realized by a convolutionless, non-Markovian quantum master equation, developed 
specifically to treat non-Markovian interactions with the environment m ■ Such an equation has been used to describe 
energy transfer in molecular systems where the electronic excitation couples to non-Markovian vibrational modes of 
the environment @1 ESI I5U- Then the dephasing contribution can be written in the form 

£ Deeha S(p(t)) = £ w + _ L ^^(t)p{t) _ ^Dt j (4) 

n 

where the non-Markovian influence of the environment is captured by time-dependent auxiliary operators A® (t ), 
which follow a separate evolution equation that is independent of the density matrix pit) (see Refs. 0, ’5Q11S])- The 
system operators 

L* = \n){n\ (5) 


couple electronic excitation of site n to the environment. We note that the coupling to the non-Markovian environment 
that is described by Eqs. Q and ([5]) not only induces dephasing, but can also induce relaxation within the vibrational 
manifold of the environment (52] . However, in the remainder of the paper we will refer to the term Eq. Q simply as 
“dephasing”, to distinguish it from the electronic relaxation (of the electronic states) that we describe next. 

Electronic dissipation/relaxation between electronic states is analogously described through 


£ Relax (p(t)) = Y, (Ll lP {t)A*\t) + A R (t)p(t)L R 1 - L*?A*{t)p{t) - p(f)H R ' \t)L* ) , 

l,n^l 


( 6 ) 


where the system coupling operators 

T R =\l)(n\ (7) 

describe relaxation from a state | n) to a state \l), and the corresponding time-dependent auxiliary operators H R (t) 
similarly include the non-Markovian effects of the coupling to the environment. 

In the Markovian limit, where the correlation time of the environment is assumed to be short compared to the 
relevant system time scales of the dynamics, the time-dependent auxiliary operators for dephasing become time- 
independent ( 501 : 

A n{t) ^\ln L n ( 8 ) 

where are the system-environment coupling parameters. For the relaxation contribution we have analogously 150] 

^ 2 (9) 

In this limit, the non-Markovian quantum master equation ([Th becomes the well-known Markovian Lindblad equa¬ 
tion l|p|. This can be easily seen by inserting Eqs. Q and ^9| in the terms for dephasing and relaxation Eqs. © 
and (roj), respectively. 
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III. PROBABILITY CURRENT ANALYSIS OF ENERGY TRANSFER 

Our goal is to analyze the energy transfer in terms of probability currents between the sites, in order to i) provide 
information about the transfer pathways, ii) reveal the respective contributions of unitary dynamics, dephasing, and 
relaxation to the transport, and iii) identify the contribution of coherence between the sites to the transport. In 
the following, we will develop this analysis based on the non-Markovian quantum master equation description of the 
previous section. 

Energy transfer between the sites occurs when the electronic excitation is transferred between the sites, driven by 
the inter-site coupling. Since the overall excitation probability in the system is conserved, "YhnPnnit) = 1, even if 
the system is open (i.e., it couples to an environment and p(t) is identified with the reduced density matrix of the 
system), and a continuity equation holds, 


^tPnnif) — ^ ( jin (1) ■ 
l^n 


( 10 ) 


where ji n {t) is the net probability current at time t that transports excitation probability from a site l to site n. If 
ji n {t) is positive, excitation is transported from site l to site n; if ji n (t) is negative, there is transport from site n 
to site l. These currents can be identified with the energy transfer between the sites. (We note that instead of the 
probability currents considered in the present work, others have considered energy currents to quantify excitation 
energy transfer [54] . However, this can be problematic in open systems, since the energy inside an open system is not 
a conserved quantity, and therefore there is no corresponding continuity equation.) 

By making use of the additivity of currents ji n (t ) between the sites, the currents between sub-complexes that consist 
of a number of sites can also be calculated. Thus, the current Jab( t) between a sub-complex A and a sub-complex B 
is given by 


Jab (t) = EE 


( 11 ) 


IgA nGB 


When the current Jab(I) is positive, there is a net flow from A to B, and when Jab( 0 is negative, there is a net flow 
from B to A. 

We now analyze the energy transfer dynamics given by the quantum master equation Q ([7]) of Section [IT] in terms 
of these probability currents. We find that this quantum master equation leads to changes of the site populations in 
time that are given by 

d t Pnn(t)= -i(n\[H,p(t)]\n) + {n \C Dephas (p(t))\ n} + (n |£ Relax (p(t))| n) 


= H in Im(pin(t)) + ]To 

l^n l^n 


E 2 Re { (1 1 (pmfj (t) - A* (t)p(t)) I n) } , (12) 


l^n 


where the elements Hi n = (l \H\ n) of the Hamiltonian are assumed to be real. Comparing Equation (12) with the 
continuity equation (10), we find that the population currents from sites l to sites n are given by 


jm{t) = J; u n nitary (i) + c pnas w + ir ax w, 


Dephas / 


(13) 


with 


•Unitary / 


h r , = 2ff in Im(pj n (t)), 

j?: phas (t) =o, 


j)n elax W = 2 He ( ^2pik(t)Af nM *(t) - A^ l lk (t)p kr 

\ k 


(14) 


(t) 


for l ^ n. The currents are a sum of the three contributions of unitary dynamics, dephasing, and electronic relaxation. 
In Equation these contributions are expressed by the coherences pi n (t ) between the sites and the populations 
p nn {t) of the sites. are the inter-site couplings and lk (t) = (l |A^(t)| k) are the site-basis matrix elements 
of the auxiliary operators describing relaxation between states n and l (* denotes complex conjugation). Giv en a 
system density matrix p(t) and the relaxation auxiliary operators A^(t) at a time t, one can then use Equation (14) 
to calculate the currents between the sites. We now discuss the three different contributions to the total current of 
Eqs. p|, p). 
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Unitary contribution: 

As can be seen in Equation (14), the unitary contribution j^ mtary (t) to the total current is caused entirely by the 
coherence between the sites, in particular by the imaginary components of this, since there are no diagonal elements 
of p(t) (site populations) present in this term. This is an important result, since it shows that without coherence 
the unitary contribution to the energy transfer would be zero. As expected, this contribution is proportional to the 
inter-site couplings Hi n . On the other hand, this unitary contribution of the currents does not explicitly depend on 
the diagonal elements of H , i.e., the energies of the sites. Nevertheless, as we will see in the following, jji mtary (t) 
does implicitly depend on the site energies through the time evolution of the coherence. The above result for unitary 
contribution J^ nitary (i) is well known, and is given for a system with nearest-neighbor interaction in, e.g. Ref. 12&j . 

An important property of this unitary contribution j^ mtary (t) is that, because the coherence pi n {t ) is limited by 
the Cauchy-Schwarz inequality 


\Pln(t )\ 2 < pu(t) p nn (t), (15) 

for the coherence and population of a density matrix p(f), the unitary contribution to the current is also bounded: 


b’jn mtary (t)| = 2|i?j r Jm(/9i n (i))| < 2\Hi n \J pu{t) p nn (t) - (Re(pi n (t))Y. 


(16) 


Thus, the larger the real component of the coherence between two sites, the more limited is the unitary current 
between these sites. It is important to emphasize the two different roles that the imaginary and the real components 
of the coherence play in the energy transfer: while the imaginary component of the coherence constitutes the actual 
transport current, the real component does nevertheless have a constraining effect on the energy transfer. Because of 
these two different roles, it can be insightful to study the imaginary and real components of the coherence involved in 


energy transfer dynamics separately. In Section III C we do this explicitly for the special case of system-environment 
coupling in the Markovian limit. 

Dephasing contribution: 


In Equation (13), the contribution of dephasing to the current is zero, regardless of the specific form of the dephasing, 
i.e., the specific behavior of the time-dependent auxiliary operators A^(t) that describe the non-Markovian influence 
of the coupling to the environment. (Here we have assumed that (n\l) = S n i, i.e., the states are orthonormal.) This 
finding is important, because it means that the dephasing does not influence the currents and hence the energy transfer 
explicitly. Thus, in a model that only includes unitary dynamics and non-Markovian dephasing, but no electronic 
relaxation between some or all of the sites - a model often applied 0 EH ES] - the currents between the respective 
sites would be given entirely by the unitary contribution j^ mtary (t). Nevertheless, from Equation (jij) it follows that 
the dephasing term £ Dephas (p(t)) does act on the coherence between the sites, via the terms 

(l I^Dephas^))! „) = (i | (p{t)A»\t) - p(t)A»\t) - A? (t)p(t) + A%(t)p{tj) | n > (1 - 8 ln ). (17) 

Thus, the dephasing term can, for example, cause the coherence to decay. Therefore, the dephasing can indirectly 
influence the unitary contribution j^ mtary (t) to the current that is driven by the coherence, and thus implicitly 
influence the energy transport. Since the action of the dephasing terms on the coherence depends on the time- 
dependent auxiliary operators A® (t ), the precise way in which coherence is influenced will depend on the details of 
the non-Markovian dynamics of the environment. 

In the Markovian limit however, the action of the dephasing terms on the coherence is simple. In this limit, where 
AnY) (Eq. (18])), the dephasing is described by Lindblad terms, 


1 I^UndhTM*))! n ) = ill" + In) Pln{t) (1 - 8 ln ) 


1 


(18) 


that simply cause the coherences between the sites to decay on a time scale given by the coupling parameters 7 ^. 
Therefore, the dephasing can diminish or inhibit the unitary contribution j^ nltary (f) that is driven by the coherence. 
It is well known that for dephasing that is fast compared to the other timescales of the dynamics, the energy transfer 
is inhibited (Quantum Zeno effect) PQ. 

Relaxation contribution: 


Equation (141 shows that the third term j/Jf lax (t), the contribution of electronic relaxation to the total current, 
depends on both coherence and populations of the sites and also on the time-dependent auxiliary operators. Therefore, 
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the contribution of this term to the total current can have a complicated dependence on the interplay of electronic 
coherence and populations with the environment. 

In the Markovian limit, on the other hand, the situation is again much simpler. Inserting A ^(t) —> ni^ni (Eq. ®) 

into the j^; elax (t) term of Equation (14), yields the Markovian relaxation contribution 


• Relax,Lindbl 

Jin 


(*) = ifnPllit) 


”/nlPnn{t)- 


(19) 


We see that this contribution now has the character of a purely classical rate equation (see, e.g., Ref. 150 ), where 7 * 
specifies the rate of electronic relaxation transport from a state 1 1 ) to a state | n ) and 7 ^ is the rate for the reverse 
process. This relaxation contribution to the excitation transfer relies entirely on the populations; coherence between 
the sites does not enter into this process. Since this current contribution caused by the relaxation terms does not 
depend on coherence between the sites, it will not be destroyed by site dephasing. Quantum master equations with non- 
Markovian dephasing terms, in combination with such Markovian Lindblad terms for electronic relaxation/dissipation 
have been used for the description of energy transfer in light-harvesting systems, where the Lindblad terms describe 
(irreversible) trapping of the energy (e.g. in a reaction center of a biological light-harvesting apparatus) [51 14) [55| [55]. 


A. Quantifying the contribution of coherence to energy transfer 


The above analysis has shown that in the absence of electronic relaxation, the excitation transfer during a time 
interval At depends entirely on the coherence between the sites, given by the off-diagonal elements of the density 
matrix p(t) in this time interval. As described by Equation (14), this constitutes the unitary contribution j^ nltary (f) 


to the current. If this coherence is zero, there will be no excitation transfer during this time interval. 

We can quantify the contribution of coherence to the energy transfer, i.e., to the currents ji n (t ), by writing the 
currents as a sum of the two contributions of the populations of the sites and the coherence between the sites 


jln(p(t)) = irW)) + jln(p nCi (*))■ 


( 20 ) 


Here the first term contains the diagonal elements of the density matrix at time t in the site basis (populations) 
and the second term contains the off-diagonal elements (coherences), i.e., p(t) = p d (t) + p nd (t). This partitioning is 
possible because the propagator V of the evolution equation Q is linear in pit). The same partitioning can be applied 
to the changes d t p nn (t ) of the populations of the sites, which are just the sums of the currents to/from the sites (see 


Eq. @) 


dtPnnit) = £ j^VW) + E^ oher (P nd W) = <« \V(p d m n) + (n\V(p nd (t))\n). (21) 

l^n l^n 


Here the contributions of population and coherence can be readily calculated by (numerically) applying the propagator 
V of the evolution equation |l]) to the diagonal and off-diagonal part of a given density matrix pit) at each point in 
time, respectively. 

In the case of non-Markovian dephasing, but Markovian electronic relaxation/dissipation, described above (see 
Eq. @), we have the simple partitioning 

C her (p nd W) = jJi nitaiy (t) and jr: P (p d (t)) = ( 22 ) 

since the unitary contribution contains only coherence and the contribution of the electronic relaxation contains 
only populations and the dephasing contribution is zero. This is an important case, since it applies to many of the 
commonly used dynamical models for energy transfer in light-harvesting systems [2 E2 E3- Again it should be 
emphasized that without coherence, the unitary current, which is often the main contribution in excitation energy 
transfer, vanishes and is therefore completely reliant on the presence of coherence. 

We note that there are other dynamical models commonly applied to describe energy transfer in which coherence 
is not explicitly present, e.g., the (generalized) Forster method or other methods using classical master equations 
(rate equations) [9] [57 ■ This is not a contradiction to the aforementioned result that coherence is essential for the 
occurrence of energy transfer (in the absence of transfer via electronic relaxation between the sites), because such 
classical rate equations can be seen as effective models derived from a general quantum master equation. When this 
reduction is done consistently, it is evident that such classical rate equations implicitly take coherence into account, 
even though no coherence terms appear explicitly in the classical rate equation. A thorough analysis of this reduction 
is given in Refs. [551 [55] . For example, a classical rate equation containing conventional Forster rates for the excitation 
energy transfer (EET) is valid only within a restricted parameter regime where the electronic inter-site coupling is weak 
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compared to the electron-vibration coupling, providing a limiting case of the general quantum master equation. Since, 
as demonstrated above, in the quantum master equation model inter-site coherence is necessary for EET regardless 
of the parameter regime, it follows that an effective Forster model implicitly takes this coherence into account, even 
though no coherence terms appear explicitly in the Forster rate equation. We note that in the literature, Forster 
transfer is often termed “incoherent” (see, e.g., Refs. mmm) since it does not contain explicit coherence terms. 
However, this is somewhat misleading, since Forster transfer does implicitly take coherence into account. 


B. Energy transfer pathways 

Aside from the analysis of the contributions of the different processes to the probability currents, the probability 
current description is also very useful to identify the pathways of energy transfer in a network of sites and to reveal 
the relative weights, i.e., importance, of the different pathways. This can provide important insight into the design- 
function relationships of biological systems such as light-harvesting complexes, inter alia. To reveal the pathways, one 
can integrate the probability currents between the sites over a certain time interval of interest, to obtain the direction 
of transport and the net amount of excitation probability that has been transported via each pathway during this 
time period. Thus, the net amount of probability APi n (At) that has been transported between a site l and a site n 
during the time At can be calculated from 


APi n (At) = 


rto+At 


dtjln(t), 


(23) 


where the currents ji n (t) are calculated beforehand through Equations (13), (141 (or through an analogous equa¬ 


tion based on a different underlying dynamical model), from a given time-dependent density matrix p(t) that was 
itself calculated beforehand through propagation with the same dynamical model, e.g. the quantum master equation 
described in Section HT1 

This approach is used in a companion paper ji] to reveal the pathways of energy transport in the Photosystem II 
super-complex, a light-harvesting apparatus driving photosynthesis in higher plants, and to gain insight into the design- 
function relationships of this important photosynthetic apparatus. To illustrate this example, one of the diagrams 
of Ref. [4] showing the transport pathways (time-integrated probability currents) between outer (antenna) and inner 
(reaction center) sites (pigments) of one of the primary photosynthetic systems in green plants (Photosystem II) 
is shown in Figure [T] The arrows in the diagram show the directions and relative importance (magnitude of the 
integrated current) of the different pathways. See Ref. (4j for a full analysis. 


C. Dynamics for Markovian system-environment coupling 


When the memory time of the environment is short relative to the characteristic time scale of the system dynamics, 
the coupling to the environment is Markovian and the energy transfer dynamics are particularly simple. Therefore, 
it is instructive to take a closer look at this specific situation. We show below that in this Markovian setting, only 
three quantities are enough to describe the dynamics. These are, i) the probability current between the sites, ii) the 
real component of the coherence between the sites, and iii) the imaginary component of the coherence. 

In the Markovian limit, the quantum master equation of Section [TT] reduces to a Lindblad equation, with Equa¬ 
tions ^ and ([9]) applying to the dephasing and relaxation terms. From Equations (13), (14), and (19), the total 
population current from site l to site n is then given by 


3% ndW (t) = -d t d ln {t) = jZ Uitary (t) + i^ aX ’ Lindbl W = Im(p Zn (t)) + - 7 * lPnn (t)) , (24) 

where di n (t ) = ( pu(t) — p nn (t))/2 is the difference of the populations of sites l and n, and V) n = Hi n is the electronic 
coupling between the sites. The unitary contribution to the population current is driven by the imaginary component 
Im (pi n (t)) of the coherence between the sites and the relaxation contribution depends on the population p nn (t ) of the 
sites. Thus, we consider the evolution equations for the imaginary and real components of the coherence separately. 
From the Lindblad dynamics (Eq. (JT|) with Eqs. (§, @ inserted in Eqs. Q, ([6])), we obtain for the imaginary 
component of the coherence 


<9 t Im (pi n (t)) = 2 Vi n di n (t) - Ai n Re(pi n (t)) - r ; „ lm(p ln (t)) + ^ {V kn Re(pi k (t)) - Vi k Re(p kn (t ))), 

k^l,n 


( 25 ) 
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Figure 1: Diagram from Ref. jl] that shows the time-integrated excitation probability currents between the pigments in a 
Photosystem II super-complex (over a time period of 1 ns), based on a numerical propagation of the electronic density matrix 
using a non-Markovian quantum master equation. (Initially only pigments 7 and 10 are excited.) The arrows show the directions 
and their thickness the relative magnitude of the integrated currents. (See Ref. [I] for details of the simulation.) 


where A i n = ( Hu — H nn ) is the energy gap between the transition energies of sites l and n. The dephasing and 
electronic relaxation rates are combined into a single decay rate T;™ = (y z D + y ° + J2k^i 7 Jk + 7For the 

real component of the coherence, we obtain 


d t Re(pi n (t)) = Ai n lm(pi n (t)) - F ; „ Re(pi n (t)) - ^ (V fe „ lm(pi k (t)) - V ik hn(p kn {t))). 

k^l,n 


(26) 


The three coupled equations (24)-rt26k completely describe the energy transfer dynamics. Since the structure of these 


equations is very simple, some qualitative observations can be readily read off from these equations. 

First, we note that the imaginary component of the coherence creates population transport via the first term in 
Equation (24), as discussed above. The population transport in turn couples back into the imaginary component 


of the coherence via the population difference di n (t) in the first term in Equation (251. The real component of the 
coherence (Equation |26|)), on the other hand, is not directly coupled to the population transport. Nevertheless, it 
influences the transport of population indirectly, because it is coupled to the imaginary component of the coherence 
via the first and last terms in Equation (26) and the second and last terms in Equation (25). These terms that couple 
the real and imaginary components of the coherence can be divided into two distinct contributions. First, terms 
depending on the local coherence, i.e., the coherence between two sites l and n that affect the population current 


between the same two sites (the second term in Equation (25) and the first term in Equation (26)). Second, terms 


that account for the non-local influence of the coherence stemming from the coupling to other sites (the last terms of 
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Eqs. (25) and (26), respectively). The first set of terms taking the local coherence into account scale with the energy 
gaps A i n between the sites, while the second set of terms describing the non-local influence of coherence scale with 
the electronic inter-site couplings V\~ n . 

We note further that the influence of the energy gaps A; ra on the population current depends on the magnitude of 
the real component of the coherence since it appears only as a factor in the R e(pi n (t)) term in Equation (25). 

Another important property of the Lindblad dynamics is that the Markovian dephasing and relaxation leads to a 
simple decay term for the imaginary component of the coherence in Equation (25), and to an analogous decay term 
for the real component of the coherence in Equation (26), with a decay time scale given in both cases by the sum T; ra 
of the dephasing and relaxation rates. 

We emphasize that when energy transfer between the two sites in a dimer is considered, i.e., the smallest possible 
coupled site system, with only two coupled sites, the non-local coherence terms are not present. Then, the reduction of 
possible coherence effects leads to significantly simpler dynamics, and for a homo-dimer, i.e., A-| 2 = 0, Equation (26) 
would even be completely decoupled from the other two equations and the real part of the coherence would simply 
decay with the decay rate given by Ti 2 . 


IV. CONCLUSIONS 

In this paper we have introduced an analysis of electronic excitation energy transport that uses excitation-probability 
currents to describe the energy transport between the electronic states. Our discussion was set in terms of a spe¬ 
cific focus on systems of coupled sites, such as molecular aggregates, light-harvesting systems, and related systems. 
We applied this analysis to an open quantum system model for energy transport, in which the dynamics are de¬ 
scribed by a non-Markovian quantum master equation that contains both terms representing unitary dynamics and 
terms that induce dephasing and electronic relaxation, which play an important role in many energy-transfer sys¬ 
tems. This probability current analysis enabled us to gain instructive insights into key features of the transport. In 
particular, it allowed us to identify and quantify the contributions of unitary dynamics, dephasing, and electronic 
relaxation/dissipation, as well as the contribution of coherence to the transport in a straightforward way. 

We found that the probability currents which represent the energy transport are given by a sum of unitary, de¬ 
phasing, and relaxation contributions. Each of these have simple forms and can be readily calculated from a given 
time-dependent electronic density matrix of the system, where this is obtained from, for example, numerical simu¬ 
lation of the quantum master equation. An important result of the probability-current analysis is that the unitary 
contribution to the current is caused entirely by coherence between the sites, i.e., the off-diagonal elements of the 
electronic system density matrix in the site basis, and not by the population of the sites (diagonal elements of the 
density matrix). This means that if there is no coherence between the sites, then the unitary contribution to the 
current will be zero. It is also noteworthy that only the imaginary component of the coherence is involved in the 
unitary current not the real component. Therefore, it is insightful to consider imaginary and real components of 
the coherence separately when studying the transport dynamics. 

Another important result is that the contribution of the dephasing term to the currents is always zero, regardless 
of specific properties of the non-Markovian environment, such as, for example, the form of the environment spectral 
density that describes the strength of the coupling between the electronic excitation and each vibrational mode of 
the environment. This means that in a model containing only unitary dynamics and non-Markovian dephasing - 
a commonly used model for energy transfer ©ESI EH - the only contribution to the overall energy transfer comes 
from the unitary contribution to the current and thus from the coherence between the sites. However, even though 
the dephasing does not explicitly contribute to the currents, it can affect the coherence that constitutes the unitary 
contribution to the currents and thereby influence the energy transport implicitly. If the dephasing is strong, it can 
destroy the inter-site coherence and thus inhibit the unitary contribution to the transport currents - this is the well- 
known Quantum Zeno effect [T] . We note that such a non-Markovian dephasing term in the quantum master equation 
can also induce relaxation within the vibrational manifold of the environment, since the environment spectral density 
can be decomposed into single vibrational modes that couple to the system and to a separate environment, leading 
to relaxation of these single modes [52] , 

We found that the contribution of the electronic relaxation to the current is more complicated and can in general 
depend explicitly on both coherence and population terms in the electronic density matrix, as well as on the non- 
Markovian dynamics of the environment. In the Markov limit, however, when the memory time of the environment is 
short compared to the other time scales of the dynamics, the electronic relaxation contribution to the current becomes 
a very simple classical rate equation term that only depends on the populations of the sites and the corresponding 
relaxation rates, but not on the coherence. We note that performing the Markov limit with respect to the electronic 
relaxation still allows the incorporation of non-Markovian dephasing - and hence of vibrational relaxation as noted 
above - caused by a non-Markovian environment [3] i52( [55] . 
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Such non-unitary electronic relaxation is often used to model radiative or non-radiative decay of electronic excitation 
or trapping of excitation, e.g. in a reaction center of a light-harvesting system mmmm- This is a separate 
process from the unitary, coherent excitation transport mechanisms between sites, which occurs, e.g., via resonant 
transition dipole-dipole interaction between molecules mmm- Thus, even in the presence of electronic relaxation, 
the actual excitation transport currents between sites often rely on the unitary contribution to the currents, since, as 
the probability current analysis shows, the third possible contribution, dephasing, does not explicitly contribute to 
the currents. In this very commonly encountered situation, the transport currents will then depend entirely on the 
coherence between the sites. Consequently, if the coherence is zero at a time t , there will be no current, i.e., no energy 
transport, at this time t and the overall energy transport is dependent on the presence of sustained coherence in the 
system. 

The finding that the energy transport may be entirely determined by the coherence is one of the key insights 
provided by this probability current approach. We emphasize that this reliance of the transport on coherence does 
not contradict the fact that the (overall) transport dynamics can often be described reasonably well by classical master 
equations, such as the rate equations deriving from Forster and modified Redfield models that contain only population 
terms and no explicit coherence terms, in contrast to the quantum master equations that contain the coherence as off- 
diagonal elements of the density matrix n mm GlZ]. This is because such classical rate equations can be interpreted 
as effective descriptions of a full quantum description, in which the effective rate equation implicitly takes coherence 
into account by virtue of derived, effective rate parameters, even though the coherence does not explicitly appear in 
the equation [22] , 

Another very useful aspect of the description of energy transport by means of probability currents is that this 
analysis can reveal the pathways of the transport and their relative importance in a large network of sites and can 
thus help to gain insight into design-function relationships of biological networks such as light-harvesting complexes. 
In this approach, the probability currents between the sites can be integrated over a certain time interval of interest to 
show the net amount of excitation probability transported via each pathway and the direction of the transport within 
this time period. The probability current analysis that we introduced in this paper is used in a companion paper [4] 
to reveal energy transport pathways in the Photosystem II super-complex that drives photosynthesis in higher plants, 
leading to new insights into the design-function relationship of this photosynthetic apparatus. 

We qualitatively analyzed the interplay between the probability current and the imaginary and real components of 
the coherence for the case of a Markovian environment where the evolution of population and coherence is described 
by a Lindblad quantum master equation. This yielded analytic understanding of the mutual influence that these 
quantities have on each other and how they are affected by dephasing, relaxation, and localization due to energy gaps 
between the sites - effects that play a crucial role in modeling and understanding excitation energy transport W- 

Finally, we point out that while the probability current analysis presented in this paper has focused on a quantum 
master equation description of energy transfer, our procedure to identify and quantify contributions to the probability 
currents that stem from the different terms of an evolution equation is more general and can be applied to a broad range 
of dynamical models and evolution equations. It can in principle be applied to any evolution equation that describes 
the time-dependence of probabilities or other conserved quantities. This includes the classical master equations that 
contain probabilities and rates, such as those derived from a generalized Forster or modified Redfield models E3H3, 
as well as other forms of quantum master equations. 
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